Casimir Interactions Between Scatterers in Metallic Carbon Nanotubes 
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We study interactions between localized scatterers on metallic carbon nanotubes by a mapping 
onto a one- dimensional Casimir problem. Backscattering of electrons between localized scattering 
potentials mediates long range forces between them. We model spatially localized scatterers by 
local and non-local potentials and treat simultaneously the effects of intravalley and intervalley 
backscattering. We find that the long range forces between scatterers exhibit the universal power 
law decay of the Casimir force in one dimension, with prefactors that control the sign and strength 
of the interaction. These prefactors are nonuniversal and depend on the symmetry and degree of 
localization of the scattering potentials. We find that local potentials inevitably lead to a coupled 
valley scattering problem, though by contrast non-local potentials lead to two decoupled single- valley 
problems in a physically realized regime. The Casimir effect due to two-valley scattering potentials 
is characterized by the appearance of spatially periodic modulations of the force. 

PACS numbers: 03.70.+k,73.63.,11.80.-m,61.72.S- 
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I. INTRODUCTION 

A single-walled carbon nanotube is a two dimensional 
graphene sheet rolled into a cylinder. The diameter of 
the nanotube is on the order of a few nanometers, and 
its length can vary from hundreds of nanometers to cen- 
timeters. Due to the small tube radius, electrons are con- 
fined in the azimuthal direction, and at sufficiently low 
energy the quantum confinement leads to an effectively 
one-dimensional electronic system. These nanotubes can 
be either metallic or semiconducting, and the low-energy 
electronic band structure can be studied using a long- 
wavelength expansion of the Hamiltonian around each of 
the degenerate Fermi points, labeled by K and K' points. 
This long-wavelength theory is given by a pair of one 
dimensional Dirac Hamiltonians. 

When a nanotube is chemically functionalized or con- 
tains defects on the tube wall, localized scattering centers 
interrupt the free motion of its low-energy charge carri- 
ers. Generally a localized defect can backscatter a prop- 
agating low-energy electron, either by large momentum 
scattering between the K and K' valleys or by small mo- 
mentum backscattering from forward to backward mov- 
ing excitations within a single valley. Superposition of 
right and left moving excitations produces various stand- 
ing wave patterns in the electron density near such a de- 
fect. 

In this paper we consider forces on the scatterers pro- 
duced by their interaction. It is easy to see that for an 
isolated scatterer, the backscattering-induced forces on 
the left and right hand side of the defect must exactly 
cancel, so there is no net force. However, for pairs of 
defects and generally for a distribution of defects at fi- 
nite density, the forces on the left and right hand sides of 
the scatterer do not balance and mediate a net force on 
each scatterer. In previous work we explored this effect 
within a single valley model for the nanotube and found 
that the scattering induced forces could be mapped to 



a Casimir-type problem, where the propagating electron 
waves provide the role of the background quantum field. 
Importantly, the spinor character of these background 
fermions admits the possibility of attractive, repulsive or 
compensated null forces on the scatterers dependingon 
the internal symmetry of their scattering potentials [l| . 

In this work we generalize these earlier results to 
study the combined effects of intravalley and interval- 
ley backscattering. This extension proves to be crucial 
for a meaningful application to the nanotube problem. 
Potentials that produce only intravalley scattering need 
to vary slowly on the scale of a lattice spacing. Yet, any 
local potential with this property degenerates to a one- 
dimensional scalar potential that cannot backscatter a 
massless Dirac particle. Thus, for a local potential our 
effect ultimately requires a significant degree of spatial lo- 
calization, and in this regime intervalley backscattering 
ultimately arises. Indeed, we find below that for local po- 
tentials there is no regime in which the force problem can 
be regarded as confined to a single valley, necessitating a 
coupled valley formulation of the scattering problem. 

By contrast, non-local scattering potentials do allow 
the possibility of only intravalley backscattering in a con- 
trolled physically realizable limit. This situation is real- 
ized most naturally for electrons coupled to slowly vary- 
ing lattice strains on a nanotube. In this paper we present 
a generalization of the formalism described in Ref. [l[ 
suitable for application to the coupled two valley prob- 
lem, and explore the forces that occur as a function of 
range and internal symmetry of the scattering potentials. 
We provide formulae that describe the electron mediated 
forces in these various geometries. Table I provides a 
compact summary of our results. 

The magnitude and sign of the interaction is dictated 
by the internal structure of the scatterers. Local poten- 
tials can describe atomically sharp impurities localized on 
a sublattice site. We find a repulsive force between lo- 
cal impurities residing on equivalent sublattice sites and 
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an attractive force between scatterers on distinct sites. 
Related results were recently shown for interactions be- 
tween impurities in two-dimensional graphene [2]. We 
also explore interactions between impurities where only 
intervalley scattering is present. Interactions between de- 
fects due to large momentum backscattering were previ- 
ously discussed in one-dimensional Fermi liquids [HOI- 
For non-local potentials we show that scattering persists 
for ranges that are larger than the lattice constant lead- 
ing to the single-valley scattering problem. The results 
we obtain for Casimir forces between non-local scatterers 
agree with our previous work. We recover the univer- 
sal distance dependent power law decay for the Casimir 
force in one-dimension. However, for local potentials, 
unlike for the single-valley problem, we also observe pe- 
riodic spatial modulations in the force due to intervalley 
scattering. 

The paper is organized in the following manner. In 
Sec. |TT] we define the geometry and derive the low-energy 
electronic structure of single- walled carbon nanotubes. 
In Sec. Mil we present scattering potentials which can 
describe impurities in nanotubes. The distinction be- 
tween relevant length scales is discussed in Sec. MI Al In 
Sec. MIBI and Sec. MI CI we discuss local and non-local 
potentials, respectively. In Sec. |IV] we outline the basic 
mechanism used to calculate Casimir forces. In Sec. II V Al 
we review our previous work of the one-valley problem, 
and show how the method is generalized to the two- valley 
problem in Sec. IIVBI Our main results are presented in 
Sec. El Casimir forces between local and non-local po- 
tentials are shown in Sec. IV Al and Sec. IV Al respectively. 
In Sec. I VII we discuss the relation of our findings to phys- 
ical adsorbates on nanotubes. The paper is concluded in 

Sec. EH 



II. SINGLE- WALLED CARBON NANOTUBES 
A. Nanotube Geometry 

In this section we describe the geometric structure of a 
carbon nanotube and introduce the notation used in this 
paper. Two-dimensional graphene is a honeycomb lattice 
with two inequivalent sublattice sites, labeled A and B, 
as illustrated in Fig. [1] There is one carbon atom residing 
on each lattice site. The primitive lattice vectors are 
ai = a(l,\/3)/2 and ai = a(— 1,\/3)/2, where a/y/3 ~ 
1.4 A is the nearest neighbor bond length. The vectors 
r's define a triad of nearest neighbor bond vectors as 
shown in Fig. [TJ 

A carbon nanotube is formed by wrapping the 
graphene sheet into a cylinder, such that two equivalent 
lattice sites arc identified. The circumferential vector 
C = n&i + ma.2, where n,m £ Z, characterizes the nan- 
otube. The xy-plane defines the lattice coordinate sys- 
tem, where bonds run parallel to the y-axis. The tube 
coordinate system is defined by Xy along the tube axis 
and x ± around the circumference. The two coordinate 



systems are related by the tube's chiral angle defined as 
the angle between x and shown in Fig. [TJ The 

coordinate transformation is given by 



cos 9 sin 9 \ / x 
- sin# cos 8 J \y 



(1) 



The circumference vectors of high-symmetry achiral nan- 
otubes that have a plane of mirror symmetry are shown in 
Fig. [TJ In armchair (6 — 0) and zigzag (9 = 7r/6) carbon 
nanotubes bonds run parallel to the tube's circumference 
and axis, respectively. 

Fixing the origin on an A site, the lattice translation 
vector Ra = + ^2&2, where ni,ri2 € Z, locates an 
A sublattice site, and the vector R_b = R4 + r locates 
a B site, where t is a vector connecting the two sublat- 
tice sites. The lattice vectors in the nanotube coordinate 
system are given by 
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where b — for i — A and b — 1 for i = B. The nearest 
neighbor bond vectors r,-'s shown in Fig. [JJ in the tube 
coordinate system are given by 

Tj = — ^z^sint'jiii — cibfl,.!' 



(3) 



where 0j = 6 - 2nj/3, and j = {0, ±1}. 

The first Brillouin zone of the honeycomb lattice is 
shown in Fig. O In graphene the conduction and valence 
bands touch at the six corner points of the Brillouin zone. 
Therefore, for undoped graphene the Fermi surface lies 
at the K and K' points. The three equivalent Fermi 
points identified by white and black circles in Fig. [2] are 
related by reciprocal lattice vectors G = mibi + m 2 b 2 . 
However, K and K' points are inequivalent since they 
cannot be connected through a reciprocal lattice vector. 
In the nanotube coordinate system the six corners of the 
Brillouin zone are given by 

4tt 
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where a = +1(— 1) for K(K')-points, and p = {0, ±1}. 
As shown in Fig. [21 the corner point K Q is a reference 
defining the chiral angle 9 between the lattice x-axis and 
the tube axis. 



B. Low-Energy Theory 

The energy band structure of graphene can be obtained 
using a tight-binding model for it electrons. Consider- 
ing nearest neighbor hopping between sites on a two- 
dimensional honeycomb lattice, the tight-binding Hamil- 
tonian for graphene is given by 



Ho = —t 



E 

Raj 



a t (R A )6(RA + 



(5) 



3 




FIG. 1: Two-dimensional honeycomb lattice with A and B sublattice sites identified. The primitive unit vectors are ai and 
&2, and t's define a triad of nearest neighbor bond vectors. A nanotube is characterized by a vector C = nai + ma.2 that 
point along the tube circumference. The chiral angle 8 is the angle between the lattice coordinate x and the tube axis x\\. The 
circumference vectors of high-symmetry achiral armchair (n,n) and zigzag (n,0) nanotubes are shown. 
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FIG. 2: The first Brillouin zone of the honeycomb lattice 
depicted relative to the tube coordinate system, where x\\ 
points along the tube axis. The six corners of the Brillouin 
zone are shown. The three equivalent K (black circles) and K' 
(white circles) points are related by reciprocal lattice vectors 
G = mibi + m2b2. The chiral angle 9 is defined as the 
angle between the tube axis and Ko in the lattice coordinate 
system. 



where t is the nearest neighbor hopping energy, and 
aJ(i)T) creates an electron on the A(B) sublattice. The 
low-energy electronic properties are found by expanding 



the tight-binding Hamiltonian around the two distinct 
^(i-T^-points to linear order in momentum k. Since the 
Fermi points are inversely proportional to the lattice con- 
stant |K| oc 1/a, the long-wavelength theory is valid for 
|k|a«l. 

The energy spectrum of a carbon nanotube is obtained 
from the graphene Hamiltonian by rotating to the tube 
coordinate system and quantizing the crystal momentum 
along the transverse direction. Single-walled carbon nan- 
otubes are either metallic or semiconducting depending 
on whether the discrete lines of crystal momentum pass 
through the Fermi points K and K' . It turns out that 1/3 
of all nanotubes are metallic, since mod(n — to, 3) = is 
a necessary condition for the six corners of the Brillouin 
zone to be allowed wave vectors. 

In our notation, the 2x2 identity and Pauli matri- 
ces {Icr,(Ji} span ^4(_B)-sublattice pseudospin space, and 
{I T , n} span the K (_RT')-point valley isospin space, where 
i = {x,y,z}. For simplicity, we introduce an operator 
which defines a rotation by an angle r\ around n in cither 
t or a space. For example, in a space this operator is 
given by 

0'(n a ,r)) = e i^n/2 0e -iA-<7r,/2_ ^ 

It is convenient to define a projection operator P~£- = 
(1 ± cr 2 )/2 which projects on a sublattice site. Likewise, 
P^- = (I T ± t z )/2 is a projection operator in the valley 
space. 

In this paper, we only consider the lowest energy band 
of metallic tubes (gapless systems) as will be explained 
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in Sec. IIIII Expanding 7i around the Brillouin zone 
corners and rotating to the tube coordinate system using 
Eq. (jTJ) , the long- wavelength Hamiltonian for the lowest 
energy band of a metallic nanotube becomes 



— iflVF 



d x 



(7) 



where Hvp — \/3at/2 ~ 0.54 eV-nm. The basis states are 
four-component spinors defining relative amplitudes at 
the A and B sites and the K and K 1 Fermi points in the 
following order (AK p , BK p , AK' p , , BK' pl ), where p and 
p' correspond to one of the three equivalent K and K 1 
points, respectively, depicted in Fig. [5] The eigenstates 



of H , f±l{x }l 



/v27r are right and left mov- 
ing plane waves multiplied by a spinor, where k is the 
momentum along the tube axis. When the chemical po- 
tential is fixed at /i — the filled Dirac sea has E = — \k\, 
and the right and left moving spinors are given by 




( ° \ 



V2 



(8) 
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functions in Eq. (fTO)) are rapidly oscillating and describe 
modulations on the scale of the atomic spacing. Since 
F(q) decreases rapidly with momentum, in the lowest 
"star" approximation [(J we keep terms in the sum of 
G„'s which connect the three K(K') Brillouin zone cor- 
ners, such that |K + G| = |K|. This approximation is 
appropriate for the range of the scattering potentials we 
study in this paper. The normalized basis functions at 
the A and B sites in the lowest "star" representation are 
given by 



v m— U,±l 



K m -r — am 

Z , 



(11) 



m=0,±l 



where z = exp(i2-7r/3). Evaluating the matrix element of 
a tight-binding potential in the lowest "star" basis given 
in Eq. (| 1 ljl and expanding to linear order in k, one ob- 
tains the low-energy Hamiltonian given in Eq. ([7]) for the 
lowest band of a metallic nanotube. 



III. SCATTERING POTENTIAL 



C. Basis States 

The eigenstates of the long- wavelength Hamiltonian in 
Eq. ([7]) are isotropic and do not depend on the crystal 
orientation of a nanotube. To include lattice anisotropic 
potentials in the theory, we reconstruct the Bloch func- 
tions from the solutions in Eq. ([7]) for the effective mass 
theory. In the kp approximation the electron wave func- 
tion near the Fermi energy is given by a Bloch function 
at the K point multiplied by an envelope function. For 
graphene, the wave function is 

*(K + k,r)= Y, e Jk - r ^ K (r)^, k , (9) 

i=A,B 

where vE^kWs are exact Bloch functions at the K point, 
and e lk ' r 0i l k's are slowly varying envelope functions @. 
Bloch states are plane waves multiplying a cell periodic 
function. Potentials which resolve the lattice structure 
couple to the lattice periodic component of the Bloch 
states. Taking the Fourier transform of the periodic part 
of the Bloch function, the sublattice basis functions at 
any of the six corner points aK p 's are given by 

= e MKp ' r Y F(\aK p + G n \)e iGn < T - Ti) , (10) 

a 

where F(q) is the Fourier transform of a localized orbital 
function, G's are reciprocal lattice vectors, and a defined 
in Eq. (|U) labels the K and K' points. The subscript % la- 
bels a sublattice site, such that ta = and Tb =t . The 



In this paper, we study Casimir interactions between 
two scatterers mediated by the conduction electrons of a 
carbon nanotube. In this section we describe the struc- 
ture of the scattering potentials used to study this prob- 
lem. We explore the dependence of the Casimir interac- 
tion on the symmetry, range, strength, and orientation 
of the two potentials. We discuss two types of potentials, 
local and non-local, which result in different scattering 
processes. 



A. Potential Range 

In our previous work we studied the one-valley scat- 
tering problem valid for potentials whose range is larger 
than the lattice constant, where intervalley scattering 
does not play a role. The 2x2 matrix structure of such a 
potential is described by its pseudospin polarization [l[. 
When the range of the potential is on the order of inter- 
atomic spacing, the two valleys are no longer decoupled 
Q. In this paper, we build upon our previous work to 
incorporate the effects of sharper potentials resulting in 
a two-valley scattering problem. When the two valleys 
are coupled, the potential is described by a 4 x 4 ma- 
trix and is characterized by both pseudospin and valley 
polarizations. 

In general, the spatial variation W of the scattering po- 
tential relevant for Casimir interactions is shorter than 
the conduction wavelength of the envelope function A, 
such that Wk <C 1. Fig. [3] shows an illustration of a scat- 
tering process. Freely propagating electrons in regions I 
and III have a wavelength A oc 1/k, and the scattering 
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region II has a width W . A potential can be described by 
delta-function as long as W -C A. The important distinc- 
tion between the one- and two- valley scattering problems 
described by the spinor structure of the Hamiltonian is 
relevant for potentials whose range is longer and shorter, 
respectively, than the interatomic separation. 



I ' * X„ 

I II III 

FIG. 3: An illustration of a scattering process. I and III define 
regions of free propagation along the tube axis. The shaded 
scattering region has a width W. A scattering potential can 
be represented by a delta-function when W is much smaller 
than A, the wavelength of the envelope function. 

We study interactions between scatterers in metallic 
nanotubes. Since Casimir interactions mediated by mas- 
sive fields are exponentially suppressed at long distances 
Q , in this paper we do not address semiconducting nan- 
otubes or scattering between bands which do not pass 
through the Fermi energy. The momentum transfer in 
the azimuthal direction between various Fermi points 
is determined by the matrix structure of the scattering 
potential V. The free degree of freedom is the longi- 
tudinal momentum, and the scattering process is truly 
one-dimensional along the tube axis. We model a delta- 
function scatterer by a one-dimensional square-barrier 
potential of the form 

V(x ll ) = V9(x ll -x 1 )0(x 2 -x ll ), (12) 

where V describes the internal structure of thepotential, 
and W — (x2 — x\) -C A is the barrier width [l|. 

In the rest of the paper, long-range potentials imply 
a range d longer than the lattice constant but shorter 
than the envelope function wavelength a < d <C A. 
Short-range potentials refer to atomically sharp scatter- 
ers whose range is comparable to or smaller than the 
lattice constant d < a. 



B. Local Potentials 



cylinder V(r) is given by 

V( X]l ,x ± ) = Ve X p{ - — sm (-^-) j t 

(14) 

where V is the potential strength, r = is the 

center of the Gaussian on the nanotube surface, R is 
the radius of the tube, and d controls the range of the 
potential. 

The matrix elements are calculated in the lowest "star" 
basis defined in Sec. Ill CI For example, the intravalley 
matrix expectation value Vaa of the potential given in 
Eq. (fT3")) evaluated in the lowest "star" basis defined in 
Eq. (fTTj) is given by 

(^(r)|^(r)|^(r)) =~ £ / d\ e -^- K ^V(v) 

m,m' 

=1 £ V(K m -K m ,), (15) 

m,m' 

where V(q) is the Fourier transform of the potential. The 
Fourier transform of the Gaussian potential in Eq. (|14| is 
normalized such that V(q) — > 1 as {q n ,q ± } — > 0. There- 
fore, V(q\\,q x ) is given by 




(16) 

where I n {x) is a modified Bessel function of the first kind. 
In the large radius limit, the Fourier transform of the 
Gaussian potential approaches the limit of a potential on 
a two-dimensional flat sheet and becomes isotropic. In 
the R^> a limit Eq. (fTB)) is given by 

V(q) = Ve-^ d2/i e- l ^°. (17) 

We define the center of the Gaussian by r„ = + 
where < v < 1, such that the potential is centered on 
either the A sublattice, the B sublattice, or along any of 
the three bonds defined by the triad of bond vectors Ti 
pointing away from r^. The total impurity Hamiltonian 
is given by 



A local potential can be represented as 

V(r,r') = V(r)S(r-r'). (13) 

We are interested in the matrix structure of the scatter- 
ing potential as a function of its range and position on 
the lattice. For example, if we consider Gaussian model 
potential V(r) = Ve~ ][r ~ r °^ d ' , then on a surface of a 



Hi = HI + HI, (18) 

where HI and H\ are 4x4 matrices containing intravalley 
and intervalley matrix elements, respectively. 

Initially, we focus on the intravalley part of the po- 
tential. Evaluating both the diagonal and off-diagonal 
matrix elements Vij 's of a local potential, the intravalley 
part of the Hi becomes 
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H% = I T ® VkP+ + VbP; 



2tt{£ - p') 



2n(p - I) 



(19) 



V(v) 



V 




FIG. 4: Potential amplitudes Va, Vb and Vab denned in 
Eq. (|20p represented by dashed, dotted and solid curves, re- 
spectively, for zero potential range d/a ~ as a function of 
v. The parameter v determines the center of the Gaussian 
potential along a bond connecting two neighboring sublattice 
sites: v = indicates a potential that is A-sublattice cen- 
tered, v = 1 yields in a B-sublattice centered potential, and 
v = 1/2 corresponds to a bond-centered potential. 

The component of the potential that points along the 
electron's propagation direction does not backscatter 
since it simply shifts the longitudinal momentum and 
can be removed by a gauge transformation [9[ . Applying 
the gauge transformation, we find that the component 
of the off-diagonal matrix elements which contributes 



backscattering is proportional to Vab sin0£, where t la- 
bels the bond where the center of the potential is posi- 
tioned. When the potential is centered in the middle of 
the bond Va = Vb, and the diagonal matrix elements 
result in a scalar potential represented by an identity 
matrix. There is no backscattering by a scalar potential 
in metallic nanotubes due to Berry's phase of the wave 
function under a spin rotation [icl |. The off-diagonal in- 
travallcy matrix elements vanish when a bond-centered 
impurity is on a bond that is parallel to the tube circum- 
ference (sin.0^ = 0). For example, in Fig. [T] the circum- 
ferential vector C labeling an armchair (n, n) tube runs 
parallel to the bonds labeled by a vector r G . Therefore, if 
the center of the Gaussian is positioned in the middle of 
any r Q bond, there will be no intravalley backscattering 
by this local impurity for an armchair tube as labeled in 
Fig. HI 

Bonds are parallel to the circumference only in arm- 
chair nanotubes, and a mirror reflection about the axis is 
accompanied by an exchange of an A and B sublattice. 
Therefore, a mirror reflection across the nanotube axis 
for armchair tubes commutes with the Hamiltonian. If a 
potential commutes with the Hamiltonian, left and right 
moving states will not mix, and there will be no backscat- 
tering. Therefore, perturbations that are symmetric with 
respect to mirror reflection about the tube axis have zero 
intravalley backscattering amplitudes pdlfr3 |. 

In the large radius R^> a limit when the Gaussian po- 
tential becomes isotropic as shown in Eq. (|17[) . the coeffi- 
cients in Eq. p^|) within the lowest "star" approximation 
are given by 



where 



2 -^ 2 d 2 /4 



V A = VI 1 + -e 



V B = vl 1 + V^ d2 / 4 



(2lIV\ f^1TV\ 

2coa {—) +cos {—J J 

/2tt(i/-1)\ /2tt(2i/+1) 
( ff(2 7 1) ) + cos(gl r ' 



(20) 



r 



|K P -K P ,| 



4tt\/3 
3a ' 



V + V (21) 



is the momentum transfer between equivalent Fermi 
points depicted in Fig. [5] in the lowest "star" approxi- 
mation. 

For short-range potentials, the matrix structure of 



7 



the scattering potential is a function of the center of 
the Gaussian potential v. A plot of the amplitudes in 
Eq. (f20|) as a function of potential center v for d/a ~ 
is shown in Fig. 2J The dashed, dotted, and solid 
curves represent V A , V b , and Vab, respectively. When 
the Gaussian potential is centered on the A sublattice 
(y = 0), there is no amplitude on the B sublattice 
(Vb = 0) and vice versa. The off-diagonal amplitude 
Vab is zero for both A (y = 0) and B (y — 1) sublattice 



centered potentials and is maximum when the potential 
is bond-centered (y = 1/2). When the potential is cen- 
tered in the middle of the bond the three amplitude are 
equal Va = Vb = Vab- For long-ranged d/a > 1 poten- 
tials, the lattice structure resolution is smeared, and Hf 
becomes a scalar potential which does not backscatter 
massless fermions. 

The intervalley matrix elements that describe scatter- 
ing between inequivalent K and K' points are given by 



J 



W? = V a t' x {z t ^ a ) ® P+ + V b t' x {z t , <j> B ) ® P- + ® (aL(2a, -4> P ab) + &S) > 



r 



(22) 



where the phases are 



K R 



A- 



K R° 



2tt/3( P + P '+£), <f AB = K-R° A -2Tr/3(p-£), andK-R^ = 
27r/3(n — m ). The intervalley scattering coefficients 



V' A ,V' B , and V' AB in the large radius limit and the lowest 
"star" approximation are given by 



J 



V 



2e -Q?rf 2 /4 



?d 2 /4 




l + 2cos(— ) 

2tt(u- 1) 



-Qt,d 2 /4 



1 + 2 cos 



1 + 2 cos 



ir(2v - 1) 



l + 2cos(— ) 



+ e -Q 2 2 d 2 /4 



r 



1 + 2 cos 



2tt(2v+ 1) 



Id 2 / 4 



1 + 2 cos 



2tt(2^- 1) 




(23) 



Within the lowest "star" there are two magnitudes of 
momentum transfer between distinct Fermi points which 
are given by 



Q 1 = |2Kp| = ^ 
3a 



87T 



(24) 



Intervalley amplitudes are equal to their corresponding 
intravalley amplitudes for atomically sharp potentials 
when d/a ~ 0. 

The intervalley amplitudes approach zero for long- 
range potentials d/a> 1, unlike the diagonal intravalley 
terms in Eq. (|20p which approach a constant. Intravalley 
and intervalley amplitudes given in Eq. ([20]) and Eq. (|23|) . 
respectively, are plotted as a function of potential range 
d/a in Fig. [5] The curves labeled Va, Vb, and V' A are am- 
plitudes of a A-sublattice centered (y = 0) potential. Due 
to three-fold rotational symmetry of the lattice V' B = 
for a ^4-sublattice centered potential. When d/a ~ 
the amplitudes for intravalley and intervalley scattering 
become equal V A = V' A — 3V, and V B = 0. The vice 



versa is true for a -B-sublattice centered (y = 1) scat- 
terer. Off-diagonal intravalley and intervalley amplitudes 
Vab and V' AB vanish for a sublattice centered potential. 
The remaining two curves are plots of off-diagonal am- 
plitudes due to a potential centered in the middle of a 
bond (y = 1/2). In general, the intervalley amplitudes 
decays slower than the intravalley ones, since Q\ < Q . 
When the potential is anisotropic for R ~ a, the relative 
magnitude of the intervalley and intravalley amplitudes 
is a function of the tube's chiral angle. 

To summarize, for a Gaussian model potential inter- 
valley scattering amplitudes decay as a function of d/a 
for all values of v and are negligible for a long-range po- 
tential. The intravalley components of a local potential 
Hamiltonian also do not contribute to scattering when 
the potential is long-ranged. When the range of the po- 
tential is on the order of interatomic spacing d/a > I, the 
potential in Eq. (fTS]) becomes a scalar and is described 
by an identity matrix I T ® I a , which produces no scat- 
tering for massless Dirac fermions [T3] ■ This holds for all 
values of v, since the position of the potential is irrele- 
vant when the potential is slowly varying on the scale of 
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V(d/a) 




FIG. 5: Intravalley and intervalley amplitudes given Eq. (|20p 
and Eq. (|23|l , respectively, due to a local Gaussian potential as 
a function of range d/a for various values of potential center 
v. The plots labeled Va,Vb, and V'a are due to a A-sublattice 
centered potential (y — 0). In this case, V B , Vab, and V'ab are 
zero. For d/a > 1, Va = Vb ~ V and — > 0. The remaining 
curves labeled Vab and V'ab are off-diagonal amplitudes due 
to a bond-centered potential (y — 1/2), which decay to zero 
for a long-range potential. 



the lattice. Therefore, only atomically sharp local poten- 
tials produce backscattering, a regime where both intra- 
and inter- valley scattering play a role. Note, within our 
model one cannot realize a local potential where only 
intravalley scattering is present. Therefore, a local po- 
tential inevitably results in a two- valley problem. 



C. Non-Local Potentials 

In this section we present an example of a one-body 
non-local potential and show that it backscatters even 
when the potential is long-ranged. We model a non-local 
potential by 



V(r, r') = v(^-) [g(r r')6(r - r' - Tj ) 



+ g(r' -r)5(r-r' + Tj) . (25) 



The prefactor V(f) depends on the average f = (r + r')/2 
of the spatial coordinates, and the remaining terms de- 
pend of the difference of r and r'. The 6- functions restrict 
the length scale of g(Tj) to the nearest neighbors. The 
quantity g(Tj) can describe, for example, local modula- 
tion of the hopping integral between neighboring sites. 
This term depends on the orientation of the j bond in a 
nanotube. 

Calculating the off-diagonal intravalley matrix element 



Vab in the lowest "star" we find 
<*»|V(r)|^(r)) = ^- £ V(K m , K m ) 



X] COS [ 



(26) 



In order to obtain the dependence of the potential on the 
orientation of the lattice with respect to the tube axis, 
we study the first three terms in the gradient expansion 
of 9(Tj) given by 



s(Tj) ~ g 



g2 



(27) 



where g is a scalar, gi is a vector, and g2 is a tensor of 
rank two. We include deviations of the hopping ampli- 
tude to zeroth order in the momentum expansion around 
the Brillouin zone corners. We fix the defect potential in 
the plane of a tube's coordinate system and obtain the 
dependence of the perturbation potential on the tube's 
chiral angle 8. 

The off-diagonal intravalley matrix elements for a non- 
local potential have terms that are non- vanishing for zero 
momentum transfer. We evaluate the m = m 1 compo- 
nent of the sum in Eq. (f2"6"| for the first three terms in 
the gradient expansion of g(r,-) shown in Eq. (j2"7| . The 
zcroth-order scalar g Q term, the average of the hopping 
amplitudes, has no off-diagonal contribution at the Bril- 
louin zone corners. The first-order term proportional to 
gi is a vector potential that shifts the electronic spec- 
trum around a Fermi point. Vector potentials that cou- 
ple to the longitudinal momentum have no effect on any 
physical properties and can be eliminated by a simple 
gauge transformation. Therefore, only the the compo- 
nent of the vector potential that shift the momentum in 
the azimuthal direction can scatter incoming states. The 
second-order term in the expansion couples to gj , a ten- 
sor of rank two. These potentials describe deformations 
such at strains, twists, and curvature. Some examples of 
such perturbations can be found in 13lfT^ |. 

Including the first three terms in the gradient expan- 
sion, the dimensionless sum over g(Tj) that enters the 
to = ml term of Eq. (|26[) is given by 



Es( t j) ; 



ay/3 _, 



+ JL=e*>[(G ±± -G l 



g± + igw 

-i(G ll± + G ±]l )} 



ge 



where we have used exp(iK r 



i) 



(28) 



and 



z±m ~ 0- The components of the two-dimensional 
vector potential gi along the tube axis and circumference 
are defined by g^ and g ± , respectively, and have dimen- 
sions of inverse length. The vector potential does not 
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depend on the chiral angle as seen in Eq. (|28|) . The com- 
ponents of the rank two tensor are defined by G^ with 
dimensions of inverse length squared. For example, the 
diagonal components G u and G ±± can result from uniax- 
ial strains along the axial and circumferential directions, 
respectively. The off-diagonal components Gu ± and G±\\ 
can represent strains such as local twists [l3lll4l | . The ten- 
sor potential preserves the symmetry of the honeycomb 
lattice since it is invariant under the transformation of 
the chiral angle 9 by 2ir/3, which is apparent in the 39 
dependence in Eq. l[2"8|). 

Gauging away the component of the potential that cou- 
ples to the longitudinal momentum, the non-local defect 
potential due to zero-momentum transfer is given by 

H 2 = VIm(g)[p-®a' y (z„,6' p )-P+®a' y (z cr ,-9 p )\. (29) 

When V^(r) is modeled by a Gaussian potential, all 
other matrix elements of a non-local potential decay 
oc exp(— Qfd 2 /4) where Qi's are defined in Eq. (|2"Tj) and 
Eq. (|24|) . Therefore, these matrix elements are paramet- 
rically smaller than the ones described in Eq. (|29p for 
non-zero d/a and will not be considered further. 

The perturbation Hamiltonian due to a non-local po- 
tential given in Eq. (|29|) is independent of the impurity 
position v and preserves the rotational symmetry of the 
lattice. The potential is non-zero for potential ranges 
that exceed the scale of the lattice. The range of this 
potential is only limited by the envelope square barrier 
defined in Sec. IIII Al Therefore, for a non-local potential 
only intravalley scattering contributes for finite range po- 
tentials, and the problem is single-valley. 



IV. FORCE CALCULATION AND 
SCATTERING MECHANISM 

In our previous work we developed a framework for 
studying Casimir forces between potentials relevant for 
the one- valley scattering in metallic carbon nanotube [l[ . 
In this paper we discuss potentials where both intra- and 
inter- valley scattering are present. In this section we re- 
view the one- valley force calculation, and then generalize 
the method to the two-valley scattering problem. 



A. One- Valley Problem 

In Ref. [l] we employ the force operator approach to 
calculate Casimir forces between one-valley scattering 
potentials mediated by one-dimensional massless Dirac 
fermions. The total Hamiltonian Ti. for the one- valley 
problem is given by 



n 



-ihv F P+ <g> a' x (z a , -9 p )d x + V{x). (30) 



The first term in Eq. (f3"0"l) is the 2x2 low-energy Hamilto- 
nian expanded around the K p point, obtained by decou- 
pling the two valleys in Eq. ([7]) . The internal structure of 



the scattering potential is dictated by its spinor polariza- 
tion. We study potentials with sharp walls and calculate 
a force as the walls becomes impenetrable. We model 
a delta-function potential by a square barrier and study 
limits of zero width and infinite potential strength. The 
potential V(x) is given by 

V(x) = Ve la ^ 2 a z e- l ' 7 ^/ 2 9(x - x 1 )9(x 2 - x), (31) 

where (j) is the spinor polarization of the potential, and 
9(x) is a step function. 

The force operator is given by 



F 



an 

dx ' 



(32) 



where x is a position. Using the Hcllmann-Feynman the- 
orem the total force is the ground state expectation value 
of the force operator summed over all occupied states. 
The force exerted on one barrier with sharp walls is given 
by 



F 



dk 



+ W/2)\ V\<S>{x + W/2)) 
W/2)\V\${x- W/2)> , (33) 



where x — (x\ + x 2 )/2 is the center of the barrier, and 
W = x 2 — x\ is its width. The wave functions in Eq. (|33|) 
are linear combinations of right- and left- moving eigen- 
states of the one-valley unperturbed Hamiltonian. The 
relative amplitudes of the propagating states are defined 
by transmission and reflection coefficients. The scat- 
tering coefficients are obtained from the transfer matrix 
&(x 2 ) = T$(xi) relating the wave functions at the two 
boundaries of a barrier. The two expectation values in 
Eq. (f3"3")) represent the difference between pressures on 
the right and left sides of the barrier. For the one scat- 
terer system the pressures on both sides of the barrier 
are equal, and the net force exerted on the scatterer is 
zero. 

A non-zero force arises from multiple reflections of 
states between two or more scatterers. A scattering pro- 
cess between two barriers due to a right-moving state is 
illustrated in Fig. [BJ The scattering potentials are la- 
beled by their spinor polarization (p. The reflection and 
transmission coefficients resulting from scattering pro- 
cesses within the same valley are shown in Fig. [BJ For 
example, Rkk the amplitude of a right-moving K state 
backscattered into a left-moving K state. 

To calculate the force between two barriers, we fix the 
position of the left barrier and differentiate the Hamilto- 
nian with respect to their separation z. The total force 
is given by 



F = 



dk, 
2tt 



2-J2\R,kk\ 2 -J2\ T * 



KK 



(34) 



where the sum is over coefficients due to right and left 
incoming states. The first term in Eq. (|34[) is an outer 
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vwv 



-w- 



fV(x) 

V 

Tj,KK 



R 



i,KK 



-w- 



T, 



KK 



> X 



Region I 1 



« — z — 

Region 



Region 



FIG. 6: A one- valley scattering illustration due to a if-point 
state incoming from the left. The two barriers of width W 
and height V are separated by distance z. Each barrier is 
characterized by its spinor polarization <j>. The scattering co- 
efficients are labeled in each region of free propagation. 



pressure pushing the barriers together, and the remaining 
terms represent an inner pressure pushing the barriers 
apart. In the two barrier system, the outer and inner 
pressures are not equal resulting in a non-zero force. 

We obtain a force whose sign and magnitude depends 
on the relative spinor polarization (50 = 2 — 4>i °f the two 
scatterers. The force between two barriers separated by 
distance z in the strong and weak strength T — VW/Hvf 
limits is given by [l[ 



F 



' 2Az 2 



1 - 3 86/tt 



r> i 



12r 2 cos((5</>)/7r 2 , r«i 



(35) 



In the r> 1 limit — -k < 86 < n beyond which the force is 
periodic. When two potentials are aligned at 56 — 27m, 
we obtain a universal attractive force for the fermionic 
Casimir effect in one-dimension. When 86 = (2n + 1)tt 
the relative spinor polarization of the two scatterers is 
antiparallel resulting in a repulsive force. The oscil- 
latory dependence on 56 persists in the weak strength 
limit. Note, it is convenient to express the force in the 
strong limit in terms of a dilogarithm function Li2{x), 
F = hv F Re[Li 2 (-e t5 ' t ')}/2Trz 2 in Eq. §5§ when T > 1. 
The results in Eq. ([35]) are plotted in Fig. \7\ 

In the r ^> 1 limit the states between the barriers are 
quantized, and the number of states changes by one when 
86 is an odd multiple of tt resulting in the cusps seen in 
Fig. [7] The weak limit does not exhibit this behavior, 
since the quasibound states between the scatterers are 
described by a continuous spectrum. 




FIG. 7: Force between two barriers as a function of their 
relative spinor polarization 8(j>. The solid and dashed lines 
represent the forces in the large and small potential strength 
limits given in Eq. (|35|l . respectively. The magnitude of the 
force in the weak potential limit, the dashed curve, is rescaled 
to r = 1/2 so the two curves can be compared. 



scattering of states between different valleys as well as 
within the same valley is present. Therefore, the poten- 
tial is described by a 4 x 4 matrix characterized by sublat- 
tice and valley degrees of freedom. The intra- and inter- 
valley matrix elements are obtained using the Bloch ba- 
sis states described in Sec. Ill CI The freely propagating 
states are eigenstates of the effective Hamiltonian given 
in Eq. ([7]). The wave functions used to calculate the force 
expectation values obtained from the Hellmann-Feynman 
theorem are linear combination of right and left moving 
states from the two K and K' points. The relative ampli- 
tudes of the propagating states are defined by scattering 
coefficients. A general expression for the wave function 



V 

-w- 



o o 



B. Two- Valley Problem 

In this section we generalize the method described in 
Sec. IIV Al to the two-valley scattering problem, where 



FIG. 8: An illustration of a scattering mechanism by a square 
barrier potential described by a matrix V and width W. 
A 4 x 4 scattering matrix is obtained by relating right and 
left moving K and K' states to their corresponding outgoing 
states. 
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in a region of free propagation is given by 



$(x) = 
+e 



OLK><t>k 



(x) =e lkx w (a K c}> 



(36) 



where </>'s are four component spinors given in Eq. I|8]). 
and a's and /3's are scattering coefficients. For simplicity 
of notation we have dropped the p and p' superscripts 
referring to one of the three equivalent corner points. The 
three K points are related by reciprocal lattice vectors, 
and physical quantities will not depend on the particular 
choice of the corner point. The dependence on p and p' 
enters only as a phase of the scattering coefficients a's 
and /3's. 

The full Hamiltonian for the one square barrier system 
is given by 



Ht =H + VQ{x l - xi)9(x 2 - x {l ), 



(37) 



where TL Q is the low-energy Hamiltonian given in Eq. ([7]) , 
V is a perturbation potential, such as 7i\ or 7i 2 described 
in Sec. lIIIi and the step functions define a square barrier. 
Integrating Eq. (|3"T|) across the barrier, the 4x4 transfer 
matrix becomes 



T = exp I - iW[P+ <g> a' x {z a , -9 p ) 

-P-^a' x (Za,e' p )][k + V}}, 



(38) 



where W = x 2 — x\ is the barrier width. 

From the transfer matrix, we calculate the scattering 
matrix. The 4x4 scattering matrix, obtained from in- 
coming and outgoing states illustrated in Fig.[8j is defined 

as 



where = (1 



a" 



t r' 
r t' 



(39) 



are right-moving incom- 



ing (i) and outgoing (o) amplitude column vectors, and 
/3's define left-moving states as shown in Fig. [5J The 
"primes" in Eq. (|3"9")1 indicate the coefficients due to the 
states incoming from the right. Each coefficient in the 
scattering matrix in Eq. (I39|) is a 2 x 2 matrix defining 
both intravalley and intervalley scattering amplitudes. 
For example, 



t= ( tKK 

\tKK> 



tfC'K 
tji'K' 



(40) 



where the diagonal(off-diagonal) terms are the intraval- 
ley(intervalley) transmission coefficients. For instance, 
tKK' is the forwardscattering amplitude of a right- 
moving K state being transmitted into a right-moving 
K' state. 

As in the one- valley problem, non-zero forces arise from 
interactions between two scatterers. An scattering pro- 
cess illustration of a left-incoming K state between two 



potentials V\ and V 2 separated by distance z along the 
tube axis is shown in Fig. [5J As before, we fix the left 
barrier and calculate the force exerted on the right bar- 
rier using the Hcllmann-Feynman theorem. The force is 
given by 



F = 



dk, 

■z— k 

2tt 



E 



Til 



(41) 



where the summations represent a sum over all reflection 
and transmission coefficients in-between the two barri- 
ers (region II in Fig. [5]) due to right and left incoming 
states, (p± k and 4>± k - Throughout this paper lower-case 
coefficients will refer to scattering by one barrier, and 
upper-case ones due to scattering by a two barrier sys- 
tem. 

The first term in Eq. (|41[) represents an outer pressure 
in Regions III of Fig. [9] due to a continuous spectrum 
of states pushing the barriers together. The second and 
third terms in Eq. (|41|) result in the inner pressure push- 
ing the barriers apart, which is obtained from the coeffi- 
cients in Region II of Fig. [9] These coefficients are given 
by 

r<=ti + ti(l-r»/ 1 )- 1 rati 



Tl = t' 2 + r 2 (l - r[r 2 ) 1 r[t' 2 
Ri=r 2 {l-r' 1 r 2 )- 1 t 1 
Rl = r[(l-r 2 r[)-H> 2 . 



(42) 



When the intervalley matrix elements are zero in one 
of the scattering potentials V, there is no forward- and 
back- scattering between inequivalent Fermi points for 
the two-barrier system. In this case Eq. (|4"Tj) reduces to 
the one- valley force given in Eq. (f34j) . 



V. RESULTS 

Using the method described in Sec. IIV1 we explore the 
dependence of the force between two scatterers on the 
matrix structure, range, and strength of the defect po- 
tentials. We distinguish interactions between local and 
non-local potentials discussed in Sec. IIIII We show that 
the Casimir force decays as 1/z 2 which is a universal 
result in one-dimension in the far field limit. However, 
we also find that in the presence of intervalley scatter- 
ing there is a spatially periodic modulation of this force. 
Our results pertain to the limit z 3> W where shape 
corrections are negligible 0] . A general solution of the 
integrals appearing in the force calculations in derived in 
the Appendix, and a summary of our results is presented 
in Table [J 



A. Forces between Local Potentials 

In this section we first consider interactions between 
local potentials. As discussed in Sec. IIII B[ backscatter- 
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FIG. 9: A right-moving state <f> k is scattered by a two barrier system separated by distance z along the tube axis. Each barrier 
has a width W, height V, and is labeled by an 4 x 4 matrix-valued potential V. Generally, each potential can produce both 
intravalley and intervalley scattering as labeled by the appropriate coefficients in each region of free propagation. 



ing from a local potential is significant for potential that 
vary on the scale of the lattice d/a < 1. Let us spe- 
cialize Eq. (|18p to describe impurities that are centered 
at either of the two sublattice sites. We first study the 
strong potential limit by fixing the area of the potential 
r = VW/ Hvf- The force is independent of the magnitude 
of the potential in the r>l limit and is relevant for the 
discussion of universal Casimir interactions. For a sub- 
lattice centered potential in the atomically sharp limit 
d/a — > intra- and inter- valley amplitudes are equal 



J 



Vi ~ V[, as shown in Fig. [5J All reflection and transmis- 
sion coefficients for such scatterers approach the same 
value in the strong potential limit, |ry| = \tij\ = 1/2 
V {i, j} — {K, K'} and are independent of the sign of the 
potential. 

Calculating the two-barriers scattering coefficients de- 
scribed in Eq. (|42p and inserting into Eq. (|4"Tj) . the force 
between two impurities centered on equivalent sublattice 
sites is given by 



Faa. 



BB 



n 



kdk 



1 - cos 4 (K ■ R ) 



1 + cos 4 (K • R ) - 2 cos 2 (K • R ) cos(2fcz) 



r 



(43) 



where R is a primitive translation vector in the tangent 
plane separating the two impurities, and z the compo- 
nent of their separation along the axial direction. The 
subscripts AA and BB imply a force between impu- 
rities which are located on equivalent sites. Applying 
Eq. (A13), the solution of the force integral in Eq. ([4"3"[) 
is given by 



k AA,BB = ~ ~-k«2 

2nz z 



cos 2 (K • R Q ) 



Unlike in the one- valley problem where the force decays 
monotonically as 1/z 2 , in addition the two- valley prob- 
lem results in a spatial modulation of the force, as ob- 



served in the argument of the dilogarithm function in 
Eq. (|44[) . The force oscillates with the period of the 
V3 x V3 superlattice indicating coupling between the 
two valley points. The force given by Eq. (f4~4")) is plotted 
in Fig. [10] as a function of z/a for an armchair tube. The 
points on the curve indicate the discrete values of the 
force in each period. The force between two equivalent 
impurities is purely repulsive, as seen in Fig. 1101 since 



(44) Li 2 



cos 2 (K • R ) > 0, where cos 2 (K ■ R ) = {1, 1/4}. 



Next, we consider interactions between impurities re- 
siding on different sublattice sites. A force between an 
A-centered (y — 0) and a B-centered (y — 1) scatterer is 
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FIG. 10: Forces between sublattice centered impurities as a function of position. The force Faa,bb between equivalent 
impurities given in Eq. (|44[1 and Fab between defects residing on different sites given in Eq. (|46[1 is plotted as a function of z/a 
for an armchair tube in the strong potential limit. The continuous limits of the force functions are shown by dashed curves in 
order to stress the periodicity of the spatial modulation of the forces. The points indicate the discrete values of the force. The 
inset shows equivalent results in the weak potential strength limit given in Eq. (|47|l . 



given by 



Fai 



hvp 

7T 



kdk 



1 - 



1 - sin 4 (K ■ R G 



l + sin 4 (K-R 



r 



2sin^(K • R Q + 6») cos(2/cz) 



(45) 



where 9 is the chiral angle of a nanotube. 
Eq. (A13) the force in Eq. (|45|) becomes 



Fab — 



2ttz 2 



Lin 



sin 2 (K • R 



Applying 



(46) 



For unlike impurities the force is purely attractive for 
all values of the chiral angle. The argument of the 
dilogarithm takes three values sin 2 (K • R D + 9) = 
{sin 2 (0),sin 2 (27r/3 + 9), sin 2 (47r/3 + 9)} which also con- 
tains v3 periodicity. The force given in Eq. (|4"6"| is plot- 
ted in Fig. [10] on a curve labeled Fab for an armchair 
tube as a function of position. 

Eq. (|46f indicates that the system is invariant under 
the rotation of the chiral angle by tt, rather than by 27r/3 
as for a defect-free lattice. This occurs because the im- 
purity is fixed on the lattice rather than on the tube's 
coordinates, and the position of the scatterer co-rotates 
with the lattice for various values of the chiral angle. 
Therefore, the three-fold symmetry in the presence of an 
atomically sharp impurity is broken. The chiral angle 



dependence appears only in the force between unlike im- 
purities, since the separation between the two defects is 
not a primitive lattice vector. The three branches in one 
period of Fab are plotted as a function of 9 in Fig. [TT] 
The figure indicates that force oscillates between and 
—irhvp/ 24z 2 for all values of K • R Q . An attractive and re- 
pulsive interaction between defects on different and same 
sublattice sites, respectively, was recently shown in two- 
dimensional graphene Q. 

Next, we study the small potential T <C 1 limit and 
compare results to the ones obtained in strong r » 1 
limit given by Eq. (|1S| and Eq. (|4*6| . We keep the first 
non-zero term in the expansion of small T and take the 
zero width limit W — > 0. The next order term in the 
small width expansion accounting for shape corrections 
is 0{W/r) [l[. For simplicity, we study the case of arm- 
chair nanotubes 9 — and find a general expression 
for a force between sublattice centered defects. The off- 
diagonal matrix elements Vab and V AB are zero for sub- 
lattice centered potentials v — {0, 1}. The force between 
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FIG. 11: The three branches in one period of Fab, a force 
between an A and a B sublattice centered impurities, given 
in Eq. (|46[1 as a function of chiral angle 9. The force is scaled 
by a factor of nhvF /2Az 2 and is found to be attractive for all 
values of 6. 



two local potentials in the f<l limit is given by 




+(r'i + r'i) ■ (v'l + r'j) cos (2K ■ r g ) , 

(47) 

where s = 1(— 1) refers to a force between potentials of 
the same(different) sign of T, and the superscripts indi- 
cate the potential describing scatterer one and two. Un- 
like in the large strength limit shown in Eqns. (|4"5|) - (|4l))) . 
the sign of the force is a function of the relative sign s of 
the two potentials in the weak limit. The sign of the force 
also depends on the relative sublattice centers of the two 
scatterers, as in the strong potential limit. Therefore, in 



the r <C 1 limit the sign of the force is controlled both by 
the sublattice position of the two defects and the relative 
sign s of their potential strength. The \/3 x V3 periodic 
oscillation persists in the small strength limit. These 
results for specific sublattice positions of the two poten- 
tials and general chiral angle are shown in Table U and are 
plotted as an inset in Fig. [TU] for an armchair tube. For 
long-range d/a > 1 potentials the force approaches zero 
for all values of v since the sublattice intravalley matrix 
elements T^'s and T^'s become equal, and intervalley 
terms T^'s and r' B 's decay to zero as shown in Fig. [5] 
This result confirms the absence of backscattering from 
an scalar potential by massless Dirac fermions. 

Although a scatterer where the two valleys are decou- 
pled cannot be realized for a local potential, a case of 
pure intervalley scattering is possible. For a local po- 
tential, when an impurity is centered in the middle of a 
bond that points along the circumference, the potential 
scatters states only between inequivalent valleys as dis- 
cussed in Sec. IIII Bl This holds because the intravalley 
part of the Hamiltonian is a scalar potential for all 
values of d/a, since Va = Vb for a bond-centered poten- 
tial, and Vab = when the perturbed bond points along 
the circumference. The intervalley amplitudes are equal 
V' A = V' B = V'ab for v = 1/2, when 9 = and £ = 0. 

In this case, the intervalley transmission coefficients 
\tKK'\ = |i_RT'_R"| — and the intravalley reflection co- 
efficients \tkk\ — Vk'K'\ = vanish. The absence 
of back- and forward- scattering within the same val- 
ley and between different valleys, respectively, by poten- 
tials that preserve mirror reflection symmetry about the 
tube axis has been also shown by Ando et al. [l5| . In 
the r 3> 1 limit, the non-zero coefficients have limits 
\rKK'\{\rK'K\) ~* 1 and \t K K\{\tK'K'\) -* 0. The phase 
of the reflection coefficients depends of the sign of T. The 
force between two potentials with only intervalley scat- 
tering contribution in the large potential strength limit 
is given by 



F e = 



2hv f 



kdk 



1 - 2 lim 



|1 + s(l - T 2) e 2i(fc2-K-R <) )|2 



(48) 



where r is the magnitude of the transmission coefficient. 
The second term in the integrand representing the inner 
pressure is fundamentally different from the ones seen 
in Eq. (|43[) and Eq. (|4"5|) . The phase that appears in 
Eq. ([48]) is associated with large momentum backscatter- 
ing. The forces shown in Eq. (|4"5)) and Eq. (|45[) involve 
two types of momentum transfer which appear as vari- 
ous terms in the equations. When both intra- and inter- 
valley play a role, there is finite transmission even in 



the strong potential limit. When only intervalley scat- 
tering is present, the strong potential limit results in an 
impenetrable wall limit since transmission coefficient ap- 
proaches zero. Therefore, the inner pressure in Eq. (|48[) 
results from resonant states between the boundaries. The 
overall prefactor in Eq. (|48p is twice the magnitude than 
in Eq. (O and Eq. ([43]) . 

Applying Eq. (A14) and evaluating the periodic part 
of the force, the solution of the integral in Eq. (|4"8"]) is 
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0- 
-0.02- 



U V ■ V V'\AAi z/a 



w s=1 



FIG. 12: Forces between local impurities where only intervalley scattering is present. The force F e given in Eq. (148 \ between 
two potentials of equal (s = 1) and unequal (s = —1) signs is plotted as a function of z/a for an armchair tube. The continuous 
limits of the force functions are shown by dashed curves. The points indicate the discrete values of the force. The inset shows 
equivalent results in the weak potential strength limit given in Eq. (I5U|I . 



given by 



-Re 



Li2(—se 



2iKR„ 



~72z~ 2 



{-3,1}, 8 = 1 

{-2,6}, s = -l 



(49) 



When only intervalley scattering amplitude is present 
the force oscillates between attractive and repulsive with 
V3 period as observed in Eq. (|49|) . The magnitude of 
the force is determined by the relative sign s of the two 
potentials. A plot of F e as a function of z/a for s = 
±1 is shown in Fig. [T3] The points in the plot indicate 
the discrete values in each period of oscillation given in 
Eq. In the small strength limit T <C 1 the force 

becomes 

F e = ^ shvp ^ 2 cos (2K • R ). (50) 

The results of Eq. (|50|) are shown as an inset in Fig. [T2j 
Although the prefactors of the force are different in the 
two limits, the oscillation between attractive and repul- 
sive persists in both weak and strong potential limits. 
Similar behavior has been observed previously in one- 
dimensional Fermi liquids where only large momentum 
backscattering is considered [KH . Refer to Table U for a 
compact summary of the main results presented here. 



B. Forces between Non-Local Potentials 

In this section we calculate Casimir forces between 
impurities described by non-local potentials given in 
Eq. ([29]) . When the range of a non-local potential is 
d/a > 0, off-diagonal intravalley matrix elements Vab 
are dominant since all other amplitudes are parametri- 
cally smaller as noted in Sec. IIII CI Therefore, a non- 
local potential can result in a one- valley scattering prob- 
lem discussed in Sec. lIV Al These potentials can describe 
modulations to the hopping amplitudes between neigh- 
boring sites. In the absence intervalley scattering, states 
are scattered only within the same K point. Therefore, 
the scattering coefficients \tkk>\ — Vk'k\ — \tKK'\ = 
\tK'K | = are zero. Likewise, the intervalley coefficients 
due to states incoming from the right vanish. Since the 
two Fermi points are decoupled the perturbation matrix 
is described by two independent 2x2 matrices in the 
sublattice er-space. 

The control parameter we vary to study interactions 
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Form 


Range 


Site 1 


Site 2 


Force (r > 1) (Eq.) 


Force (T < 1) (Eq.) 


T ^^.ol /TT/i 11 Qh 

ijocal (l^q. lioj) 


d/a ~ 


v = 
y = 1 


!/ = 
1/ = 1 


|^Li 2 [cos 2 (K-R )] @U 


- £ S f cos 2 (K-R ) (S3 


K = 
1/ = 1 


V = 1 

v = 


27rz j ^2 [— Sin ^JV ■ rt ■+- t>)\ [I40D 


sin (K • Ho + 6*) (|47[) 






r, "f r> r [ r „• c , p 2 » ki m1 (HaTi 


shv F T 2 ^/rii/ -r> \ JF7TI 

-^2 — cos(/±v ■ rt ) l|oUI) 


(9 = Q,£ = Q) 


d/a £ 1 


Any 


Any 








Non-local (Eq. [29j 


d/a / 


\V\>1 
\V\<1 


|V]>1 

|y]<i 


_2*2f d52]) 


(EH 


\V\>1 
\V\<1 


|y]<i 

|V|>1 




(35} 



TABLE I: A summary of results described in Sec. IV Al and Sec. IV Bl The first group present results of forces between local 
potentials. The remaining rows show results for forces between non-local potentials, where the dependence of the force of the 
relative sign s of the potential strength \V\ is stressed. 



between two non-local defects is the sign of the potential 
V. We assume that the dimensionless quantities Im(<?)'s 
defined in Eq. (|28j) are equal for the two barriers. In the 
strong potential T 1 limit the magnitude of the non- 
zero scattering coefficients approach |ri<-_K-|(|rK''j(-/|) - * 1 



J 



and \tKK\{\tK'K' |) — > 0. We calculate the interaction 
between two barriers with the same and different signs 
of r = Vlm(g) / hv f ■ Applying the one-valley force re- 
sult given in Eq. (j34|) . the force between two non-local 
potentials becomes 



Fo 



2hv f 



kdk 



1 - 2 lim 



o |l + s(l-r 2 )e 2ifez | 2 



r 



(51) 



where s is the relative sign of the two potentials. The 
integrands in Eq. 1|48|) and Eq. (fSTjl are equivalent except 
for the phase exp(2iK • R Q ) appearing in Eq. (|48p . This 
phase associated with large momentum backscattering is 
absent in Eq. (|51[) since there is no intervalley scattering 
present by potentials given in Eq. (f2"Tj|) . 

The solution to the integral in Eq. (|5"Tj) is shown in the 
Appendix. Applying Eq. (A14) the force is given by 



hv F irhv F 
F 2 = — jii2 ~s = — — 



-1, 



1 



(52) 



2, s = -l 



The result in Eq. (p)2")) shows that there is an attractive 
force between two scatterers with equal sign of T (s = 1) 
and a repulsive force between defects of unequal sign of 
r (s = —1). The relative sign of V is analogous to the 
difference between the spinor polarizations 5<p of the two 
scatterers discussed in Sec. lIV Al Potentials of equal sign 
(s = 1) refer to the case of parallel scatterers 8<fr = 0. 
Two potentials of opposite sign (s = —1), on the other 
hand, refer to the case of anti-parallel scatterers 6<f> = it. 
The results in Eq. ([52")) are consistent with the force in 



the T > 1 limit of Eq. ([35]) [l|. Likewise, F 2 in the r«l 
limit agrees with Eq. ([55]) . The magnitude of the force is 
larger than the result in Eq. ([55]) by a factor of two since 
we are including fermions from the two K(K') branches 
of carbon nanotubes. These results are shown in TableHJ 
Intervalley scattering becomes important for non-local 
potentials when Im(g) = for d/a < 1. A few example of 
such defects are a vector potential with a zero component 
along the tube axis (g^ — 0), a tensor potential for arm- 
chair tubes and zero twist (9 = and G ]]X = G ±ll = 0), 
or a tensor potential for zigzag tube with zero uniaxial 
strain (9 = n/6 and G tn = G±± = 0). The effect of in- 
tervalley scattering on the Casimir force is discussed in 
Sec. IV Al in the context of local potentials, and the same 
physics apply for the case of non-local potentials. 



VI. DISCUSSION 

Defects or impurities on a carbon nanotube can 
backscatter electrons either through intravalley or inter- 
valley scattering processes. In general both channels are 
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present with their relative strengths determined by the 
range and symmetry of the scattering potentials. The 
models we present here provide a framework for un- 
derstanding the backscattering-induced forces on these 
species. The signature of intervalley scattering is a spa- 
tial modulation of the scattering- induced forces. By con- 
trast intravalley scattering mediates a force that can be 
either attractive or repulsive, but has a strength that 
decays monotonically as a function of increasing separa- 
tion. Interestingly, in all cases where the interaction is 
described by a local potential, the scattering problem is 
inevitably multivalley in character, and the energy and 
force of the species oscillate as a function of separation. 

The long-range interaction between multiple scatterers 
might lead to complex phase structures. It was suggested 
by Shytov et al. 2] that interaction between adatoms ab- 
sorbed on the graphene lattice can result in defect aggre- 
gation and inhomogeneities on the lattice. 

The scatterers we describe in this paper can be phys- 
ically realized by various atomic and molecular species 
adsorbed on the tube wall. These range from cova- 
lently bound atoms and molecules [l6lfl7| , to more weakly 
bound metallic species [13]. The range of the scattering 
potential is determined by the size of the absorbed species 
relative to the lattice constant. The symmetry of the po- 
tential is determined by the spatial variation of on-site 
energies and by the modulation in the intersite hopping 
amplitudes produced by these species. 

Covalently bound species provide the most natural 
candidates for the strongly coupled local potential models 
described in section ITlI Bl Here, the on-site potential bar- 
rier at an adsorbed site can be as large as 5 eV enforcing 
an effectively hard wall boundary condition on the elec- 
tronic wave functions. In this regime the results of section 
IV Al can be used to provide a bound on the electron- 
induced force. For example, the maximum attractive 
force between two scatterers in the impenetrable wall 
limit leads to an interaction energy of E c — —Trhvp/12z. 
With Hvf ~ 5.4eV- A for nanotubes this gives an en- 
ergy of 2.8 meV at a range z = 50 nm. Note that its 
spatial form follows the same scaling law as the Coulomb 
interaction between uncompensated charges, but it is re- 
duced by a factor irhvp/l^e 2 ~ .1. Thus, for charge 
neutral dipoles p — es whose electrostatic interactions 
scale as E d ~ —p 2 /z 3 — — (e 2 /z) x (s/z) 2 , they are dom- 
inated by the Casimir interaction in the far field z > 5s. 
Similarly, this one-dimensional Casimir interaction com- 
pletely dominates the familiar van der Waals interactions 
between charge neutral species that are mediated by the 
fluctuations of the exterior three dimensional electromag- 
netic fields. 

The weak coupling limit is relevant to the interactions 
of less strongly bound species, such as metal atoms or 
molecules bound by ir stacking interactions, e.g. ben- 
zene. Here the energy scale for the local potential is more 
modest, of order 1 eV which, assuming a range of order 
a graphite lattice constant, corresponds to a dimension- 
less coupling parameter T ~ 0.5. In this weak coupling 



limit Ei = —Hvf^ 2 /2irz a local potential of V ~ 1 eV 
results in 0.4 meV at a distance of z — 50 nm. Though 
weaker, this interaction still decays slowly as a function 
of distance (ex 1/z) and will also dominate the electro- 
static interaction between charge neutral dipoles in the 
far field. 

In this weak coupling regime, strain induced couplings, 
represented by non-local scatterers can be comparable in 
size. Assuming a linear scaling of intersite hopping am- 
plitudes with bond lengths following dt/dt ~ 4eV/A a 
bond length change of 0.2 A and a potential range on the 
order of the lattice constant, this gives a dimensionless 
potential strength of T ~ 0.37 and a weak coupling inter- 
action E n i — — Hvf F 2 /irz, we find 0.2 meV. These are of 
the same order as the forces produced by local potentials 
in the weak coupling limit. 

For adsorbate-induced potentials, it is difficult to re- 
alize a regime where the scattering is dominated by po- 
tentials with solely a nonlocal form. Thus, one concludes 
that intervalley scattering and a residual spatial oscilla- 
tion of the force is a generic property of inter adsorbate 
interactions mediated by the propagating electrons. It 
may be possible to quench the intervalley channel by ap- 
plication of a magnetic field along the tube axis which 
would have the effect of introducing a gap at either the 
K or K' point and isolating the effects of intravalley scat- 
tering. We also note that strains can be engineered into 
these structures by application of mechanical stresses, 
and this might provide an avenue for realizing the pre- 
dictions of the nonlocal model. 



VII. CONCLUSION 

In this paper we show that interactions between scat- 
terers in metallic carbon nanotubes results in a one- 
dimensional Casimir problem. We generalized our previ- 
ous work which includes the one-valley problem of nan- 
otubes, to incorporate the effects of intervalley scatter- 
ing. We show that local potentials in nanotubes produce 
a two-valley scattering problem. The decoupling of the 
two valleys is not possible for a local potential since the 
range must be atomically sharp in order to produce fi- 
nite backscattering. Local potentials whose spatial ex- 
tent is beyond the lattice constant result in scalar poten- 
tials which do not backscatter massless Dirac fermions. 
Non-local potentials, on the other hand, can result in a 
decoupled valley scattering problem. Intervalley scatter- 
ing amplitudes are parametrically smaller for finite range 
non-local potentials. Therefore, we formulate a physi- 
cally realizable potential which reduces to the one- valley 
scattering problem. 

We study forces between two scatterers mediated by 
the propagating electrons of metallic carbon nanotubes. 
For interactions between both local and non-local po- 
tentials we find a universal 1/z 2 power law decay for a 
one-dimensional Casimir force. However, for local poten- 
tials, where intervalley scattering plays a role, we also 
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observe a position dependent periodic modulation of the 
force. The signs and magnitudes of the forces are not 
universal and are controlled by the internal symmetry of 
the scattering potentials. 
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APPENDIX A: FORCE INTEGRALS 

In this appendix we provide a derivation for the inte- 
grals that appear in the calculations of Casimir forces for 
r> 1. Although, a cutoff function is introduced in order 
to control divergences appearing in the integral, we show 



that the final result is cutoff independent. The class of 
integrals found in this paper have a general form 



1 f 00 
F = — kdk 

n Jo 



I-P 2 



1 + p 2 ± 2pcos(2kz + ip) 



(Al) 



where z is the impurity separation along the tube axis. 
The integrand in Eq. (|A1|) can be represented in terms 
of a Poisson kernel 



i-p 2 



+ p 2 + 2ps cos(g + ip) ' 



(A2) 



where s = ±1, and q = 2k z. Introducing an exponential 
cutoff function, the integral in Eq. (|A1|) becomes 



F = lim 



qe 



1-P sp (q,<p) dq, (A3) 



Since the Poisson kernel is 2tt periodic in q, the integral 
can be expressed as an infinite sum times an integral over 
a region of [0, 2tt]. Rewriting Eq. (|A3[) we obtain 



F = lim 



™ 4vrz 2 J [ 



1 - Psp{q,<fi) 



dq 



^(q + 2nw)t 



— (j,(q-\-2rnr) 



(A4) 



Expressing the sum in terms of a geometric series and 
separating terms constant in q, the series in Eq. (|A4[) to 
O(fi) is given by 



- p{q+2nir) 



2tt 



2tt 



(l-e- 27r ^)2 \-e-2-Kn 



dfj, \ 1 - e- 2n ^ 

g (27T-g) 
47T 



(A5) 



The first two terms on the RHS of Eq. (|A5[) diverge in 
the limit fj, — ► 0, but vanish when integrated over q since 



1 - P sp (q,(p) dq = 0. 



(A6) 
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To verify that the above statement is true in the case of 
p — > 1 we express the Poisson kernel in terms of a delta 
function 

oo | 8(q - q n ), s = l 
limP sp fe^) = 27r^^ (A7) 

% ~ ?n)> S=-l 



n=0 



where q n = n(2n + 1) — ip and = 27m — ip. Insert- 
ing Eq. (|A7p into Eq. (|A6|) , we find that there is either 
one 5-function in the range of integration [0, 2tt] or two 5- 
functions at the two limits of integration, each contribut- 
ing half the area. Therefore, in both cases the integral 
over the series of ^-functions yields a factor of 2tt, which 
is consistent with the result in Eq. (|A6[) . Note, in the 
p — ► 1 limit Eq. (|A3j) can be solved using a generalized 
Abel-Plana formula which provides a finite expression for 
a difference between an infinite integral and an infinite 
sum fl9l ). 

Combining the above results and noting that the third 
term in Eq. (|A5[) is cutoff independent, Eq. (|A4[) becomes 



F 



1 



q(2n-q) l-P sp (q,<p) dq. (A8) 



We use the following identity to solve the integral in 
Eq. CM): 



(A9) 



7T- / f{x)g{x)dx = V f{n)g{-n), 

Z7T /n — 



where the "hat" indicates the Fourier series of the original 
function. The Fourier series of the Poisson kernel is given 
by 



e in{q+ip) p\n\ 



s = -l 

P sp {q,v)= 

n=-co I (_\\n e in{;q+tp) p\n\ ^ s= \ 

. lA ' m 

The Fourier series of the other term in Eq. (|A8|) is given 
by 



27T v — v ( 



n— — oo 



Using the results from Eqns. (|A"9"] ) - (|Alip . Eq. 

comes 



F 



2ttz- 



■E 

n=l 



cos(mp) p n 

( — l) n cos(rnp)p n 



s = -1 
s = 1 



(All) 



1 be- 



(A12) 



Eq. (|A12[) is a general result which can be applied to all 
the integrals encountered in this paper. The series above 
can be represented in terms of dilogarithm functions. For 
example, 



and, 




(A13) 



cos(n^) 



s = -1 



Re[Lt 2 (-ae**')] = ]T 

n=l 

where Lii(x) is a dilogarithm function 



n=l | cos(rty) „ _ 1 

5 



(A14) 



In Sec. IV Al we calculate forces between two local sub- 
lattice centered impurities. The solution of Eq. (|43|) for 
interaction between defects residing on equivalent sites 
is Eq. (|A13j) . where ip = 0, with p = cos 2 (K ■ R Q ) and 
,s = —1. The result for the force integral in Eq. (|45[) . ap- 
plicable to interactions between impurities centered on 
inequivalent sites, is Eq. (IA13|) with p = sin 2 (K • R D ) 
and s = 1. 



The integral in Eq. (jAlj) can also be related to in- 
tegral in Eq. (|48|) for a force between two local poten- 
tials where only intervalley scattering plays a role, and 
Eq. (|5Tj) for interactions between non-local potentials. 
The limit of zero transmission r — > is equivalent to 
p —> 1 in Eq. (|A"Tj) . where p = \/l - t 2 . Writing Eq. ([48]) 
and Eq. (|5ip in a general form in terms of p we obtain 



F = — 



1-2 lim 



1-P 2 



p^l |1 + S/9 2 eI (2fez+ V )|2 

p'-^I 1 + p 4 + 2ps cos(2fcz + tp) 



1 — lim 



(A15) 



where we have ignored the prefactors. The right-hand 
side of Eq. (|A15|) is equivalent to Eq. (|A1|) in the limit 
p — * 1. Therefore, the solution of Eq. (|48|) is given by 
Eq. (|AT4| for </3 = -2K • R . The solut ion to Eq. (HU) is 
obtained by setting <p = in Eq. (IA14|) . 



